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In this paper we study decoherence in the quantum walk on the line. We generalize the method 
of decoherent coin quantum walk, introduced by Brun et al [Phys. Rev. A 67, 32304 (2003)]. Our 
analytical expressions are applicable for all kinds of decoherence. As an example of the coin-position 
decoherence, we study the broken line quantum walk and compare our results with the numerical 
one. We also show that our analytical results reduce to the Brun formalism when only the coin is 
subjected to decoherence. 
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I. INTRODUCTION 

The quantum walk (QW) is the quantum analogue of 
the classical random walk (CRW). Notable differences 
between the QW and the CRW are the quadratic depen- 
dency of variance on the number of steps and the complex 
oscillatory probability distribution in the QW instead of 
the linear variance dependency and the binomial proba- 
bility distribution in the CRW. These differences between 
the QW and the CRW have been used to present sev- 
eral quantum algorithms in order to solve some specific 
problems [J, 0, S, S, |, 1 with performances better than 
the best known classical versions. Recently Childs has 
shown that the universal computation can be performed 
by QW and it is, therefore, another witness of QW 
importance. Two types of QW have been introduced as 
the quantum mechanical extension of CRW: discrete Q 
and continuous time @. Both have the same result in 
our problems but finding the relation between them was 
an open problem for several years. Recently this rela- 
tion has been found [IS EH- Another classification of 
QW is based on the network over which the walk takes 
place. The QW on a line is the simplest possible configu- 
ration d, El, ES E3 , but other to pol ogies such as cycles 
US El, two-dimensional lattices [HI, El, El, SS SI, or 
n-dimensional hypercubes [H, [23[ have also been inves- 
tigated. 

Quantum entanglement is another important property 
of quantum mechanics which has recently attracted much 
attention in view of its connection with the QW. The 
effect of entanglement of the coin subspace on the QW 
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[HI EH, l24l ] , en tan glem ent between the coin and position 
subspaces l2a l26j and the QW as the entanglement 
generators [23, (28| are examples of these studies. Beside 
of all theoretical studies, experimental implementation 
and realization of the QW is another interesting subject 
for researchers [H, US M, 11, H, H . 

In the experimental implementation, the environment 
effects will be so important, because in practice, the 
preparation of pure quantum states without interaction 
with the environment is impossible, and the environment 
can disturb the quantum states and fades the quantum 
properties. Therefore it is very important to formulate 
and quantify the influence of decoherence on the QW and 
several valuable researches have been done about the de- 
coherent qw us u us is H3, n, is m. Due to 

the complexity of the analytical calculations, in most of 
the previous studies the numerical calculations have been 
used in investigation of the effects of decoherence on the 
QW. Although analytical expressions have been driven 
for some particular cases such as the QW with a coin 
subject to decoherence [13] or the weak noises case fill ], 
the general analytical formulas with the wide range of 
use are not found yet. The aim of this paper is to gener- 
alize the approach of Brun et al [37| and show that such 
generalization is applicable for all kinds of decoherence 
in the one-dimensional QW. 

This work is organized as follows. Section II gives a 
brief review on the one-dimensional QW and decoher- 
ence. Section III is devoted to drive analytical expres- 
sions for the first and second moments in the presence 
of general decoherence. Our results have then been used 
in section IV to analyze the coin-position decoherence for 
which separation between the coin noise and the position 
noise is impossible and, as an example, an analytical cal- 
culation on the broken line noise has been presented. An 
analysis of the coin decoherence for which only the coin is 
subjected to decoherence is also presented in this section 
and it is shown that our results reduce to the Brun et 
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al formula for the coin dccohcrence. We summarize our 
results and present our conclusions in section V. 



II. BACKGROUND 

From the various methods of studying the effects of 
dccohcrence on the quantum systems, the Kraus repre- 
sentation is one of the widely accepted method [42|. Let 
us define Hw as the Hilbert space of our system (Walker) 
and He as the Hilbert space of environment spanned by 
{|e„)}™ =0 where m is the dimension of the environment's 
Hilbert space. In practice, it is not possible to completely 
isolate the system from the environment, therefore in or- 
der to find the time evolution of the system we should 
consider the time evolution of the whole system (sys- 
tem+environment) , and obtain the state of the system 
by tracing out over the environment's degrees of freedom, 
i.e. 



p sys = Tl'Env (Uptf) 



(1) 



Here U acts both on the system and environment Hilbert 
spaces. Without loss of generality we assume that the 
state of the whole system is p = po (£> \envo}{envo\. So we 
can write Eq. JTJ) as 

m m 

Psys = E ( e »l U \ env o) Po (env \ C/ f |e„) = 2J E n p El 

ra=0 n=0 

(2) 

where E n = (e n \U \envo) , n = 0, 1, • • • ,m, are the so 
called Kraus operators. These operators preserve the 

m 

trace condition, i.e. E^En = I. 

n=0 



A. Decoherent one-dimensional QW 

QW is the quantum version of the CRW where instead 
of the coin filliping, we use the coin operator to make su- 
perposition on the coin space, and instead of the walking 
in CRW we use the translation operator to move quan- 
tum particle according to the coin's degrees of freedom. 

In one-dimensional QW we have two degrees of free- 
dom in the coin space H C7 spanned by {\L}, \R}}, and infi- 
nite degrees of freedom in the position space H pi spanned 
by {\i) i = —00, • ■ ■ , 00}. The whole Hilbert space of the 
walker is defined as the tensor product of the coin space 
H c and the position space H p , i.e. H = H p ® H c . In 
one step of the QW we first make superposition on the 
coin space with the coin operator U c and after that we 
move the particle according to the coin state with the 
translation operator S as follows 



U w = S.{I®U C ) 



(3) 



where 



S = E \ x + ( x l ® \ R ) ( R \ + \ x ~ 1 ) ( x \ ® \ L ) ( L \- ( 4 ) 



Therefore the quantum walking is defined by 

\H>(t+l)) = U w \H>(t)) - |*(t)) = C/,„ t |*(0)). (5) 

In ideal case that one can isolate the system from the 
environment perfectly, the time evolution of the system 
takes place coherently and the state of the walker re- 
mains pure after t steps of the walking. But in the real 
world where the interaction between the system and the 
environment is unavoidable, the purity of the system is 
decreased and the state of the walker becomes mixed. In 
this case we should consider the evolution of the whole 
system, i.e, system+environmcnt, defined on the Hilbert 
space 



(G) 



By definition of the Kraus operators, one step of the walk- 
ing can be written as follows 



p(*+l) 



J2E nP (t)El 



(7) 



For t steps, we can write 

m mm 

p(t) = E - E E E nt ...E n2 E ni p(0)ElEl 2 ...El. 

Tit — 712—0 71-1— 

(8) 

It is worth to note that the Eq. (jSj) is general and the 
Kraus operators E n include the whole information of all 
types of evolution. It follows therefore that the coin op- 
erator, the translation operator and the environment ef- 
fects all are embedded in E n and we didn't assume any 
restriction yet. 

The Kraus operators satisfy an important constraint 
known as completeness relation [42| which arise from the 
fact that the trace of p (t + 1) must be equal to one, i.e. 



I = Tr (t + I)) = Tr E nP {t)E n ^j 
= Tr (^EjE^t^j 



(9) 



Since this is true for all p then 

J2EjE n = I. 



(10) 



To make any progress we should therefore find the Krusc 
operators for our system defined in Eq. and use 

the Eq. |(5J) in order to obtain the final state p(t). Ev- 
idently, the E n are operators that act on the system 
(coin+position) Hilbert space, and therefore we can write 
the general form of E n as follows 



x.x' i,j 

= EEE4?Lk+o (*i®ii> 01 



(ii) 



x I i,j 
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where x, I = — oo, •■•,00 and i, j — {L, R}. 

In the following section we show that a reasonable sug- 
gestion on the coefficient aS n ' of Eq. (fTTj) enables us to 
derive useful analytical expression for the first and second 
moments of position. 



B. Environment effects (decoherence) 

In this section we briefly discuss about the interpreta- 
tion of noisy evolution and the Kraus representation. As- 
sume that the pi is the probability that the ith unknown 
reason affects the state of the system. In the other word, 
the evolution operators A, act on the system with the 
corresponding probabilities Pi . 

Now we are interested in the walker, though the uni- 
tary evolution takes place in the Hilbert space of the 
system+environment, given in Eq. ([6]). Furthermore the 
environment could be the rest of the universe, so deter- 
mining the exact reasons needs the investigation of the 
whole universe which is impossible. 

Fortunately, in practice, understanding of the exact 
form of the environment is not required. Let us assume 
that we have r different operators Ai, i = 1, • • • , r where 
each of them acts on the system with the probability Pi . 
If the environment be in the state |ej) then the operator 
Ai acts on the system. Therefore we can imagine the r- 
dimensional Hilbert space for the environment, spanned 
by {|e;)}r=i! and the following initial state for the envi- 
ronment 

|env ) = VPi |ei) + y/p~2 |e 2 ) + ••• + ^Jp~r\e r ) ■ (12) 

It is clear that the probability of finding the environ- 
ment in the state je,) is pi. Therefore we can write 
the unitary transformation of the whole system (envi- 
ronment+system) as follows 

U = |ei) (e 1 \®A 1 + \e 2 ) (e 2 \®A 2 +...+\e r ) (e r \®A r . (13) 

So the Kraus operators will be 

Ei = ( ei \ U \env ) = jpiA,. (14) 

Accordingly, the Eq. ([7]) gives the density matrix after 
the first step as 

r r 

P ' = j2 E if> E \ = Y.p* A *p A l ( 15 ) 

i=l i=l 

As we see, this is exactly our expectation from the 
influence of noise because it gives p' as a mixture of 
different evolutions with the corresponding probabilities 
Pi. With Eqs. (TT2"]) and (fT5|) we can represent the effects 
of noise by the Kraus operators. 



III. FORMALISM 

The Fourier transformation is a powerful tool for an- 
alytical investigation of the one-dimensional QW, where 
is defined as follows 

/dk 
-e-^|fc). (16) 

Using this transformation we can write Eq. (1111) in k 
space as 

& =EE^,y// _^ e -- e --( fe - fc ') |ib) (k'\ « \i) 0'|. 

(17) 

(n) 

In the following we assume that the coefficients a x [ i ■ are 
not dependent on the coordinate x. This means that the 
probability of translation in the position space depends 
only on the length I of the translation, not on the posi- 
tion x where the translation takes place. In view of this 
constraint on the , we are able to derive an analytical 
expression for the first and second moments of the posi- 
tion, which is applicable for a large family of decoherence 
in the one-dimensional QW. In this regime, the Eq. (jTTJ) 
takes the following form 

^ =EE a S / j (k - k') \k) (k'\ ® \i) 01 

l i,j 

(18) 

where we have used the orthonormalization relation 

e- ix ( k - k ') = 2tvS (k - k') . (19) 

X 

By integration on k' and changing the order of integration 
and summation we get 

/dk 
— \k)(k\®C n {k) (20) 

where 

I i,j 

Now by writing the general form of po in the k space as 

Po = JJ I*) (k'\ ® M (V'ol (22) 

the first step of walking from Eq. (J7J) becomes 

p' = //^i*x* , i®E c '«(*)i^)^i c ^(*') 

f f tikiik' 

= JJ i^rl fc H fe, l® A.fc'IV'o) (Vol (23) 

where Ck,k' is a superoperator defined by 

U.m>6 = E C n (k) OCl (k') . (24) 

n 
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Therefore after t steps, the state of the walker and the 
probability of finding the walker in position x are, respec- 
tively 



and 



P(t) = 



p(x,t) 



dkdk' 
4tt 2 



|fc)(fc'|®4 >fe ,|^ )(Vo| (25) 



dkdk' 
4tt 2 

dkdk 1 
4tt 2 



(x\k}(k'\x)Tr(£l k ,p ) (26) 
e~ ix ( fe '- fe )Tr(4 fc ,p ). 



Note that the completeness relation given in Eq. (|TCT 
implies that 



£c£(fc)<?»(fc) = J. 



(27) 



Using the orthonormalization relation (TT9|> . we get for 
the moments (x) and (x ) 



dfcdfc 



,d 2 5(k' -k) 



Tr(4 fc 'Po) 



In order to carry out these integrations we need the fol- 
lowing equations 



^(c it ,6)=T,(£^6cl(^ (33, 



It follows from this condition on the coin operator that 
the superoperators Ck,k are trace preserving, i.e. 



Tr (£ fcjfc 0) = Tr f ]T C„ (fc) OC+ (fc) 



(28) 



Tr^Ct(fc)C n (fc)6J =Tr(o 



and therefore 



Tr 



(4,0) = Tr (6) 



(29) 



which are true for arbitrary operator O. 

Evidently if we know the operators C n (k), the state 
and the probability distribution of the system can be ob- 
tained by direct using of Eqs. (|2"B"|) and (|2"6")l . But unfortu- 
nately it is these equations which arc difficult to handle, 
and therefore we turn our attention on the moments of 
this distribution. The mth moment of the probability 
distribution p(x,t) is defined by 



x m p(x,t). 



(30) 



Inserting Eq. (26]) in Eq. (J30j) we get 

r^$> m 1 1 dkdk' e-* ( fc '- fc )Tr(4 fc ,p ). 



4tt 2 



(31) 



where according to Eq. ([21 



dC n (fc) 
dk 

dCj {k') 
dk' 



= -iEE^Sj^wo-i (34) 



Since the superoperator Ck,k' ac ts on the density matrix 
which is positive and Hermitian, we can write the Eq. 
(l33l) as follows 



— Tr (£fefc'P) = Tr (Qkk'p) (35) 
^■Tr (Ckk'P) = Tr (Qh'kp) 



g^O^^^OCKk'). (36) 



Finally by carrying out the integrations of Eq. ([52]) and 
putting everything together, we arrive at the following 
relations for the first and second moments 



where 
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(•A 



w dk 



m — 1 

„tt „ t m — 1 

/ ^ E E Trj^/:™- 1 '- 1 (Wt'" 1 |Vo)(^o|) +e fc £r m ' _1 (^r'" 1 IV'oXVol) } 

w m — 1 m'= 1 



(37) 



w dk 
.„ 2tt 



m— 1 



where we have defined 



dGl 



k.k' 



dk 



k'=k 



^ dC n (k) ~ dCUk') 
dk dk' 



(38) 



k'=k 



and have used Cj~ = £k,k and Gk = Gk,k for the sake of 
simplicity. It is worth to note that thus obtained mo- 
ments in Eq. (|37[) are a generalization of the moments 
given by Brun et.al in [371 ] . in the sense that they are 
applicable for all kinds of decoherence. 




n-1 n n+1 



IV. CALCULATION AND RESULT 



(c) 



In this section our task is to show that the method is 
applicable for all possible kinds of decoherence. To this 
aim, in the first subsection we consider the broken line 
decoherence as an example of the coin-position decoher- 
ence for which we cannot separate the coin noise from 
the position noise. In the second subsection, we restrict 
ourselves to the case that the coin is subjected to deco- 
herence but the position is free of noise, and show that 
our method reduces to the result of Brun ct al given in 




A. coin-position decohrence 

In this subsection we investigate the effects of coin- 
position decoherence on the QW. As an example of our 
approach we find the diffusion coefficient of the one- 
dimensional QW in the presence of broken line noise. 
This model has been studied by Romanelli et al [4fj| . 
where the authors have calculated the diffusion coeffi- 
cient numerically They have introduced broken line de- 
coherence for the one-dimensional QW in such a way 
that the connection between neighbors of current posi- 
tion breaks with probability p, i.e. the walker doesn't 
walk with probability p. It turns out that with proba- 
bility (1 — p) 2 both links are connected and the walker 
proceeds normally. On the other hand with probability 
p(l — p) (res. p 2 ) one line (res. both lines) is broken and 
the walker is returned back (res. is stopped) (see FigfT]). 



FIG. 1: (Figure from [4(j) Possible situations for site n of the 
line when there are: (a) no broken links, (b) the link to the 
left of the site is broken, and (c) both links are broken. The 
arrows indicate the direction of the probability flux associated 
to the L,R components. 



Let us rename S in Eq. (fj| to S\ as the translation 
operator for the case that we don't have any broken line. 
We also define the translation operators S2 , 63 , S4 for 
the cases that the left, right and both lines are broken, 
respectively, i.e. 

Si = E k + 1) ( x \ ® \ R ) ( R \ + elSl \ x - 1) (*l ® \ L ) ( L \ 

S 2 =J2\x + l) (x\ ® \R) (R\ + e 1 ^ \x) (x\ ® \R) (L\ 

53 = E \x) (x\ ® \L) (R\ + \x - I) (x\ ® \L) (L\ 

5 4 = t\x) {x\®(\R) (L\+e^\L) (R\). 

X 

(39) 

The situations corresponding to the above translation op- 
erators are demonstrated in FigfT] Note that we have 
used the general phase for each translation operator. 
Later we will see that we need these generic phases in 
order to have the correct physical phenomena. 

As mentioned before, the environment state deter- 
mines the broken line situation. For instance if the envi- 
ronment is in the state |ei), then S\ acts on the system 
and so on. Therefore, according to the Eq. (TT^J) we can 
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define the initial state of the environment as follows 



\env ) = {1-p) |ei) + s/pjl - p) (|e 2 ) + |e 3 )) + p|e 4 ) . 

(40) 

Also from Eqs. ([3J and (fT5|) , the general unitary trans- 
formation of the system-environment is as follows 

U = |ei) (ei| ® 5i (/ (8) fl) + |e 2 ) (e 2 | <g> S 2 (/ ® if) (41) 
+ |e 3 ) (e 3 | ® S 3 (I ® H) + |e 4 ) (e 4 | ® &t (I ® i?) 

where we have assumed the Hadamard walk, i.e. U c = H. 
Therefore according to Eq. (fT3"]) = Si (I ® if), and we 
find from Eq. (jT4j) the Kraus operators as 



Ei = (l-p)^2\x + l){x\<g)\R){R\H 



E 3 = 



E x 



+e Wl \x - 1) (x\ ® \L) (L\ H (42) 

+e^ 2 |a!) (ar| <g> (i| (43) 

+e ?;e3 |<b-1) (x| (g) |L) (L\H (44) 
p£ |x) (x| ® (L| + e"« |L) F)(45) 



All of these operators are in the form given by Eq. (fTTj) . 
For instance, the coefficients a*- 1 -* for .Ei are 



(l) _ (i-p) 

l,R,R ~ V2 
(1) _ (1-P) 



„(1) _ -jg, (1-P) 

'1) _ pig] 

~ 6 Ta- 



rt 



where can be used in Eq. (|2"Tj) and obtain Ci as 



/A- 



_ (1-p) / e 



- ik 



V e n.*+oi) _ e i(k+e 
Similarly, we can find the other Ci as follows 

-ik 



c 2 = 
c 3 = 



(1-p) 







g-ifc _ e i 





(48) 



(1-p) ^ 

1 + e *( fc + e 3) 1 - e i ( k + 3) 



c 4 = -4 



i -i 



These C, must satisfy the completeness relation given 
by Eq. (|2"T|) . After some calculations we find that for 
any 6x,6 4 but only when 2 - 3 = tt, the Eq. Ip7|) 
is satisfied. In the following, for simplicity, we choose 
0i = 6> 3 = 6*4 = and 2 = tt. 

According to Eq. (|24[) . with these Ci in hands, we 
can find the superoperator £k,k'- But as we know, we 
have to calculate the mth power of this superoperator 
which is not an easy task to handle. To get around this 
complexity, we follow the method of Ref. [37} and use 



the affine map approach. Since the Ck is linear we can 
represent it as a matrix acting on the space of two-by-two 
operators. To do this we note that one can represent any 
two-by-two matrix by a four-dimensional column vector 
as 



O = r Q I + ri<7i + ?'2cr 2 + r 3 a 3 




(49) 



where 



1 



Tr a, Q 



(50) 



Here we have defined cto = /, and Oi (i = 1, 2, 3) are the 
usual Pauli matrices. 

Now, in order to find the superoperators £ k ,G k , Q\ and 
J7it , we represent the action of them on an arbitrary two- 
by-two matrix O as follows 



c k o 



/I 





\ 


fr 








e / + p 2 


n 








-f + P 2 e 




Vo 


1 - 2p 


-2.g -2h / 





(51) 



Q k O = 




i(p-l) 


i(J 




fr 





f 




n 





e 


l) 







-h 




\ r 3 



(52) 



(46) 



(47) and 



J k O = 



1-P 





\ 


fro 










ri 








-r-i] 


C 2 





1-P 


0/ 


\ r 3 






Qio. 





(53) 



(54) 



where the last is obtained, simply, from the Hermiticity 
of the Pauli matrices. Here e, /, g and h are functions of 
p and k, defined by 



e(p, k) 
f(p,k) 

g(p, k) 
Hp, k) 



(p- l) 2 sin(2fc) 
(p- l) 2 cos(2/c) 
p (1 — p) sin (k) 
p (1 — p) cos (fc) . 



(55) 



With this representation we can calculate the moments 

given by Eq. ([37]). Since Tr \C k O) = 2r and the C k is 

trace preserving, so it doesn't change ro- Therefore the 
only nontrivial result arise from the following three-by- 
three submatrix 




(56) 
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For finding (x) in Eq. (|3T[) we should take trace, accord- 
ingly, only the first row of Gk is important 



where a 2 



Clearly the time independent 



{X t 



-1 

7T 



dfc((p-l) .9 /i) 



E M r _1 



m — l 




We can use the geometric progression to simplify the 
summation in this equation. 



term doesn't contribute to D, therefore we focus on find- 
ing the second moment in Eq. (|37p . To begin with, let 

(x 2 ) t , i.e. 



f ^EM^r^oWoi)}. (60) 

J -* /7F m=l 



-1 

7T 











9 h) 


\l-Mu)- 1 ' 




Since only the first row of J7fc (Eq. ([55])) makes 








result in the trace, therefore 






(58) 





Note that we omit the Ml in this equation because all 
eigenvalues of Mk obey < |A| < 1 and M| — > in long 
time limit. All t dependence has vanished, therefore in 
long time limit the first moment is time independent. 

The matrix (/ — Mk) is exactly invertible and calcula- 
tion of (x)t is straightforward but our interest is finding 
the diffusion coefficient with below definition 



•>-* m =l 




(l-p)t 



(61) 



D 



1 , 9(7 2 

- hm — — - 

2 t^oo at 



(59) 



where we put rg — 1/2 because of the normalization con- 
dition of l^o)- Now we turn our attention to the first 
term of (x 2 )t in Eq. (|37l ). i.e. 



dk 



t m—l 



K E E Trj^xr™'" 1 (GkCf- 1 ivox^oi) +a fc £r m ' _1 O&t' -1 immi)} 

^ m=lm'=l 

/ l/ 2 

/ *EE( d H)j»)£! 

* m=lm'=l 



■»m— m — 1 
-fc 




r 



(62) 



where in the second line we use the facts that only the Therefore we can easily write 
first row of Gk and G\ makes nonzero result in the trace 
and that the first row of Gk is pure imaginary. 



Since the Ck is trace preserving and leaves tq un- 
changed, we write 







/l/2\ 


















\ r' 3 J 



(63) 



Also from Eqs. (|52|) and (|54|) we find the exact form of 
Gk - Qk as 



G k -0t)O = 2i i Q Q Q 








(p-l) 



Gk Gt) C^'- 1 





f f(k,p,P) 


\ 


(!)■■ 
















V p-i 


) 



(65) 



Now putting everything together, we get for the second 



moment (x 



(x 2 ) t = (l-p)t + 



dk 



T(k,p,t) (66) 



where (we drop out arguments for simplicity) 

/ 



t m— 1 



r=E E (° (p- 1 ) 3 h ) c 



m—m'—l 
k 



m—l rn' — l 



o 



Vp-i 



(67) 
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Note that the 



/ (k,p, r'j doesn't appear in T, since 
/ f(k,p,P) \ 



o (p-i) g h)cr m '- x 






V / 



= 0. (68) 



Same as before we can use the matrix M k for calculat- 
ing Eq. (Ell) 



t m—1 

E E m^'- 1 

_m=lm' = l 




r=((p-l) g h) 



where 

t m— 1 

E E M r ra " 

m— 1 m' — l 



By inserting this into Eq. (|6"5|) , the calculation of (x )t 
will be straightforward. 

Our interest is finding D here, so only the time de- 
pendent terms will be important for us. By inserting 
the time dependent term of Eq. (|70]) into Eq. ((66)) and 
definition of D given in Eq. (|59|) . we have 

1 , d° 2 
D — — hm -—— 

2 t^oo dt 



(I - M k y l E {i - MJT) (70) 

m—1 

(I - MkY 1 [t - (I - Mu)- 1 M} . 



(l-p) + 



(71) 



The matrix (/ — Mu) 1 is exactly invertible but it's ele- 
ments are a little long to appear here. Fortunately this 
matrix appear in the form of 



FIT 

£L = ((p-l) fl h){I-M k )- 1 







p-1 



(72) 



i.e., it is sandwiched between two vectors. It is not diffi- 
cult to calculate Eq. (fTTj) to have 



0.4 



0.3 



0.2 0.4 0.6 0.8 1 



= * (p) 



(73) 



where 



and 



i^(p) = -{l-(l-p)I(x)(l-p)} (74) 



7(x) 



cos (k) + x 



-7T 2tt cos (fc) x + cos (k) x + 2x 2 — 2x + 1 



FIG. 2: The K(p) versus p 



The expression of D (p) in Eq . (|73|) is exactly same as the 
numerical result of Ref. [40| in which the authors have 
proposed a constant value for K and have estimated K m 
0.4 by linear regression. But our analytical calculations 
show that the coefficient K is a function of p (Eq. ([74]) ) 
and ranges from 0.19 to 0.5 as p goes from to 1 (see Fig. 
[5]). The numeric prediction of Romanelli et al implies 
that for p=0.44 the diffusion coefficient is exactly the 
same as the unbiased classical random walk, i.e. 1/2, but 
our analytical calculations indicate that this probability 
is less than the numeric prediction and is about p=0.417. 
This means that for p < 0.417 the diffusion coefficient is 
greater than the classical one and the quantum walker 
spreads more faster than the classical one. The origin of 
this difference is behind of the fact that the numerical 
method can never be applied in problems with infinite 
steps. 



B. coin decoherence 

In this subsection we show that in the coin decoher- 
ence, the general formalism of Eq. ([571) reduces to the 
formalism of Bran ct al presented in 



37] 



Let us suppose that before each step of walking, the 

(n) 

operators D c acts on the coin subspace with probability 
p n - Then according to Eq. (fT5|) the A n are 

A n = S.{I®H).(l® =S.(l® HD^j . (76) 

By plugging it into Eq. (|14[) and using the explicit form 
of 5" we get 

En = V^E \ X + l ) fcl ® 1^) ^1 HD c n) 

X 

+ \x-l){x\®\L){L\HD l f\ (77) 
(75) Now by writing the operator HD^' in the basis 
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{\R),\m 



HDM= Y, ^\r)(s\ 

r,s={R,L} 



(78) 



the Eq. (|77p takes the following form 



E n = VKj2 E t£ |aJ + l)<s|<8|ie>H 

+ 7 g |a:-l><z|®|L>< S |.(79) 

This equation is exactly in the form of Eq. (fTTj) with the 
specified values for Z, i.e. 1 = 1,— 1, hence 



» 



» a (») 



» 



=7«;i . aii,L, s =7i,,i' (80) 
Therefore according to Eq. (|2"Tj) we have 

G»(fc) = E TSe-^I^OI+T&^I^OI (81) 



and 

dC» (fc) 
dk 



-i Y ^y ik \R){j\-^ k \L){3\ 

j={L,R} 

-iZC n (fc) . (82) 



And finally according to the definition of in Eq. 



Q hJl fO = -iZ[C k . w O 



(83) 



Clearly ^. fe ,0 = /(£/. /.■ O ) Z, and the J k in Eq. 
will be 



J k O = 



dk 



= Z[C k O)Z. (84) 



k'=k 



Inserting Eq. (|84|) into the last term of Eq. ([37]) we get 



^ Z7F m=l 



77 dfc 
ff dfc 



^Trlz^^rr 1 !^)^!)]^} 



m—l 
t 



7T 2tt 

_7r m=l 



X;Tr{^|^)(Vo|} = t 



(85) 



Finally, putting Eqs. (|83[) and (|85|) in the moment expres- 
sions of Eq. (|37p . reduces them to the same expression as 
introduced by Brun ct al [37T ] for only coin decoherence, 



(86) 



^EE Tr{z£r m '(^A n 'l^)(^o|)}+Tr{z£r m '((A n 'l^)(^|)^)}. 



m=l m' = l 



r 



Therefor the first and second moments given in Eq. ([37]) 
are general in the sense that thay can be applied for 
any types of decoherence, including coin-position deco- 
herence, coin decoherence and position decoherence. 



V. SUMMARY AND CONCLUSION 

We have present an analytical method for calculat- 
ing the first and second moments of the one-dimensional 
QW, which are applicable for all kinds of decoherence. 
Our method is a generalization of the Brun et al method, 
introduced in [371 ] . where there the authors have pre- 



sented the analytical expressions for the first and second 
moments in the presence of coin decoherence and have 
shown that the transition from quantum to classic hap- 
pens even for the weak coin noise. 

We have made exact calculations on the broken line 
decoherence as an example of the non-separable decoher- 
ence of the coin-position system. Our analytical calcula- 
tions are in consistent with the previous numerical result 
[ZoT ]. We have shown that the K in the diffusion coeffi- 
cient D is, indeed, a function of p, and ranges from 0.19 
to 0.5 as p varies from to 1. 

Furthermore our calculation shows that for p = 0.417 
the value of D is 1/2, exactly the same as the unbiased 
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classical random walk. This specific value of p is critical 
in the sense that the transition from quantum to classic 
happens at this point. In other word, the quantum walk 
with non zero broken line probability, spreads faster than 
the classical random walk if the probability of broken 
line is not higher than 0.417. For higher values of p, the 
lines are broken too frequently and this will prevent the 
full diffusion. The critical p, estimated by the numerical 
calculation of Ref. is a little higher than the one 

obtained here analytically. The difference arises from the 
non-ability of numerical method in simulation of the QW 
with infinite steps. 

We have also studied the case that only the coin is 
subjected to decoherence and have shown that the ex- 
pressions for the first and second moments lead to the 
Brun et al results presented in [37| . 

The long time behavior of moments and, of course, the 
variance can be used as a qualitative measure of classi- 
cality. Therefore these calculated moments provides a 
powerful tool for investigating the one-dimensional QW 
in the presence of any kinds of decoherence. For instance, 
finding the particular cases for which dependency of the 
variance on the time remains quadratic, i.e. the decoher- 
ence does not affect the quantumness, is an application 
of these moments. Brun et al [33| have shown that the 



quadratic time dependency of the variance is preserved 
for multiple-coin QW subjected to the coin decoherence, 
but their calculation is restricted to the coin decoherence 
only. In this paper this restriction is relaxed and the 
model allows the presence of any kind of decoherence. 

These formulas are useful also in the experimental im- 
plementation of the quantum walk. The application of 
our results for two kinds of decoherence, i.e. the coin- 
only and the coin-position noises, have been done in this 
paper. The application on the position decoherence have 
been also tested but the full calculations were so long to 
appear here. The results of our calculations on the tun- 
nelling effect, for which the position is subjected to deco- 
herence, and it's consistency with the previous numerical 
work [3l| confirm the usefulness of these expressions for 
the position decoherence too. Other aspects of the tun- 
nelling effects are also under consideration. The paper, 
therefore, can be regarded as a further development in 
the study of decoherent one-dimensional QW. 
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